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_ Abstract. We study the coherent temperature and polarization patterns produced 

in homogeneous but anisotropic cosmological models. We show results for all Bianchi 
Q^, types with a Friedman- Robertson- Walker limit (i.e. Types I, V, VIIq, VII/^ and IX) 

to illustrate the range of possible behaviour. We discuss the role of spatial curvature, 
^ ■ shear and rotation in the geodesic equations for each model and establish some basic 

c/^ I results concerning the symmetries of the patterns produced. We also give examples of 

the time-evolution of these patterns in terms of the Stokes parameters /, Q and U. 
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1. Introduction 



Precise measurements of the temperature anisotropics of the Cosmic Microwave 
Background (CMB), particularly those from the Wilkinson Microwave Anisotropy 
Probe (WMAP) [H[2j, form the sturdiest foundations of the current ("concordance") 
cosmological model describing a universe dominated by cold dark matter (CDM) and 
a cosmological constant [3], and therefore known as ACDM for short. It is an essential 
component of this model that the primordial metric perturbations that gave rise to the 
galaxies and large-scale structure we observe around us today should be Gaussian and 
statistically homogeneous (i.e. stationary) fUEllG] and this means that the temperature 
fluctuations observed in the CMB should be Gaussian and statistically isotropic. 

Detailed analysis of the WMAP data has shown that any departures from the 
standard framework are small and of uncertain statistical signiflcance. Although 
some anomalous behaviour has been reported [7J, there remains no clear evidence 
of primordial non-Gaussianity, but there are several indications of departures from 
statistical anistropy across the CMB sky. Among the interesting phenomena revealed by 
detailed analyses of the pattern of CMB temperature fluctuations are an extremely cold 
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spot P, |9[ Uni im 1121 IISI E], unusual alignments between large-scale harmonic modes 
of the temperature pattern [Ia[l6l[17l[l8l[l9l[20l[2ll[2a (sometimes dubbed 

"The Axis of Evil" ) and a global hemispherical power asymmetry [26l EH ESI ESI [30j . 

One must not get carried away with the these features because they are - almost 
without exception - based on a posteriori evaluations of statistical significance. If the 
occurrence of an anomalous feature ^4 in a given model M is P{A\M) 0.001 (say) then 
that does not necessarily mean that the probability of the model given that the anomaly 
is observed P{M\A) is similarly low. One has to take into account all the other data that 
are not anomalous before deciding on a true measure of the significance of a departure. 
Attempts to do this rigorously have generally confirmed that current observations are 
not sufficiently compelling to suggest that the standard model needs to be abandoned 
[SH [32l [33l [34] . Moreover, other analyses suggest the further possibility that the WMAP 
temperature ffuctuations may be affected by foregrounds or other systematic problems 
[35l [36l [371 [38] . Even slight effects of this type could seriously hamper our attempts to 
uncover evidence of physics beyond the standard model. 

But, although the evidence for the examples of global asymmetry discussed above 
is by no means overwhelming, taken together they do at least suggest the possibility 
that we may live in universe which is described by a background cosmology that is 
globally anisotropic, i.e. one not described by a Friedmann-Robertson- Walker (FRW) 
model and it is therefore incumbent upon us to consider alternatives to the standard 
model in order to learn best how to use the data to confirm or rule out variations on 
the standard cosmology. 

The approach we follow here is to study the Bianchi models, i.e. cosmological 
models based on exact solutions to Einstein's General Theory of Relativity involving 
homogeneous but not necessarily isotropic spatial sections. The Bianchi classification 
groups all possible spatially homogeneous but anisotropic relativistic cosmological 
models into types depending on the symmetry properties of their spatial hypersurfaces 
[39l[40]; we discuss these models in more detail later in the paper. 

The Bianchi models are not particularly strongly motivated from the point of 
view of fundamental physics, but do nevertheless represent a promising and potentially 
profitable first step away from the standard cosmological framework. For example, 
it has been known for some time that interesting localized features in the CMB 
temperature pattern can occur in Bianchi models with negative spatial curvature 
[4n [42l [43l [44l [45l [46]. The physical origin of such features lies in the focussing effect of 
space on the geodesies that squeezes the pattern of the anisotropic radiation field into 
a small region of the sky. The observed lack of large-scale asymmetries in the CMB 
temperature has in the past been used to place constraints on the global rotation and 
shear allowed in Bianchi models [Ml [Ej, [48]. More recently, however, attention has 
shifted to the possibility of using the additional parameters available in such models 
to reproduce a cold spot such as that claimed to exist in the WMAP data. Since we 
now know that our Universe is close to isotropic, attention has focussed on the subset 
of the Bianchi types that contain the FRW model as a limiting case. The model which 
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appears to best able to reproduce the anomalous cold spot is the Bianchi VII/^ case 
|49l [501 ISTl [521 [53 ] , although Blanch V also has negatively curved spatial sections and 
can therefore, in principle, also produce localised features [55]. 

However, as well as forming distinctive features in the temperature pattern, 
anisotropic cosmological models also generate characteristic signatures in the polarized 
component of the background radiation. Thomson scattering generates polarization 
as long as there is a quadrupole anisotropy in the temperature field of the radiation 
incident upon the scattering particle. In the concordance cosmology the temperature 
and polarization patterns are (correlated) stochastic fields arising from their common 
source in scalar and tensor perturbations arising from infiation. In a Bianchi cosmology, 
however, the patterns are coherent and have a deterministic relationship to one another 
owing to their common geometric origin. It has recently been shown [55l [56l [57] that 
only in special cases are the properties of the polarization field produced in Bianchi 
Yllh consistent with the latest available WMAP polarization data [58] because such 
models generally involve a large odd-parity (B-mode) [59[ [60] contribution that exceeds 
the experimental upper limit. For a discussion of CMB polarization in the context of 
the so-called Axis of Evil, see ref. [6T] . 

But their ability to produce localized features is not the only reason to be 
interested in the temperature and polarization patterns produced in Bianchi models. 
For example, history provides a host of connections between these models and the 
understanding of the interaction between electromagnetic radiation and gravitational 
waves [621 [631 [611 [Sa [661 [67]. 

Our aim in this paper is to give a relatively gentle introduction to the Bianchi 
models and show how to calculate the temperature and polarization properties of the 
radiation field expected to arise within them, using methods outlined in a previous short 
paper [55] . Our aim in doing this is neither to provide an exhaustive set of alternatives 
to the standard cosmology nor to perform a detailed statistical analysis of the patterns 
we calculate. Instead we plan to elucidate some of the general properties of the radiation 
field in anisotropic cosmological models belonging to the different Bianchi types. 

The outline is as follows. In the next section we describe the Bianchi classification 
in general terms for the benefit of the non-expert. In Section 3 we introduce the specific 
formalism of Bianchi cosmological models, i.e. exact solutions of the Einstein equations 
involving space-times with symmetries described by the various Bianchi types. In 
Section 4 we explain how we compute the radiation patterns in these models and then, 
in Section 5, show examples of the results obtained. We summarize and present our 
conclusions in Section 6. 

2. Bianchi Spaces 

In the Friedman models on which the standard Big Bang cosmology is based, 
hypersurfaces of constant time are defined to be those on which the matter density is the 
same throughout space. We can construct a more general definition of homogeneity by 
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requiring that aU comoving observers see essentiaUy the same version of cosmic history. 
In mathematical terms this means that there must be some symmetry that relates 
what the Universe looks like as seen by observer A to what is seen in a coordinate 
system centred on any other observer B. The possibly space-times consistent with this 
requirement possess symmetries can be classified into the various Bianchi types, which 
we now discuss. For more details, see [68l l69] . 

The Bianchi classification is based on the construction of space-like hypersurfaces 
upon which it is possible to define at least three independent vector fields, that 
satisfy the Killing's equation 

6;i+6;i = 0. (1) 

The vectors that satisfy this are called Killing vectors; the semicolons denote covariant 
derivatives. Suppose the Killing vectors are denoted where Greek indices can run 
from to 3. The commutators of the are defined by 

[^a, = - CfdCa = C^f^^s, (2) 

where the are called structure constants. These are antisymmetric, in the sense 
that, 

One can understand how the structure constants arise by considering symmetry 
transformations. In three dimensions, spatial homogeneity is relates to the existence of 
three independent sets of curves with tangent vectors An infinitesimal symmetry 
transformation takes an arbitrary point P (with coordinates Xi) to the point with 
coordinates Xi + where 

5x = i^5t (4) 

for (^p some linear combination of the defined at P. The same symmetry transformation 
takes the tip of an arbitrary infinitesimal vector ('p(^5 at P to a new position ^Y'+Qj^Ss^t 
at P' . This means that 

Cv^s Cp'^s = (pSs + ^p+Cp6sSt - ^p5t. (5) 

It is now possible to compare the transformed vector (p, with (^p/, the 'actual' vector 
field defined at P^ The difference between the two vectors is the Lie derivative of ( 
taken along ^ : 

C^C = lim = e • VC - C • Ve (6) 

^ st^o 5t 

This tells us, roughly speaking, by how much we must turn a vector after it is carried by 
a symmetry transformation from P to P' in order for it to point in the same direction 
as it did before the transformation. 

Instead of choosing an arbitrary vector we can now take one of the ^c^, and instead 
of an arbitrary direction we transform it along another of the The type of space is 
specified by the Lie derivatives obtained for this type of operation: 

A/^Ca = C^f^^s] (7) 
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Table 1. The Bianchi types shown in terms of whether the various parameters used 
to construct the classification are zero, positive or negative; the designation of Class A 
or Class B depends on whether a = 0, or not, respectively. The parameter h is defined 
hy h = o? /nin2' The spaces of particular interest in this paper are I, V, VIIq, Vll/i 
and IX because these contain the isotropic FRW spaces as limiting cases; the others 
are shown just for completeness. 



since it is defined as the difference between two vectors the Lie derivative must itself be 
a vector and so it can be expressed as a hnear combination of the three independent 
basis vectors. This gives us the structure constants we described above. 

The set of n Killing vectors will have some n-dimensional group structure, say G^, 
that depends on the equivalence classes of the structure constants C^^. This is used to 
classify all spatially homogeneous cosmological models. The most useful form of this 
classification proceeds as follows. On any particular spacelike hypersurface, the Killing 
vector basis can be chosen so that the structure constants can be decomposed as 

= ^aM^^"^ + ^}^oc - S^a^, (8) 

where Sa/B^y is the total antisymmetric tensor and 5f is the Kronecker delta. The tensor 
n^^ is diagonal with entries, say, ni, n2, and 77.3. The vector aa = (a, 0,0) for some 
constant a. All the parameters a and can be normalized to be ±1 or zero. If 
an2ns = then 712 and can be set to ±1 and a is then conventionally taken to be 
^J\h\ where /i is a parameter used in the classification. The possible combinations of rii 
and a then fix the different Bianchi types according to the table. 

It is interesting also to think about the generality of the different types. This can be 
expressed in terms of the dimension p of the Bianchi group, which gives the dimension 
of the orbit of (7^^ as a subset of all 9 of the distinct components. The Killing vectors 
must satisfy the Jacobi identities, so 

=0. (9) 

This implies that 

n^^a/s = 0, (10) 
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Table 2. Dependence of the scalar curvature R on the parameters a and rii involved 
in the construction of the Bianchi classification for those models with an FRW limit. 



SO that the orbits of any particular group type are at most six- dimensional. The 
(isotropic) spaces that feature in the Friedman models have Gq symmetry groups with 
Gs subgroups, so that the zero curvature {K = 0) FRW model can be thought of as 
a special case of Bianchi Types I or VIIq. Likewise the open {K < 0) FRW model is 
a special case of types V or Vllh- The closed FRW case {K > 0) is a special case of 
Bianchi type IX. 

We are interested in cosmological models that are close to the completely isotropic 
case described by the FRW metric, but not all the Bianchi types contain this as a special 
limiting case (40J|39j. Those that do are types I, V, VIIq, VII/^ and IX; we do not discuss 
the other cases any further in this paper. Bianchi I and Bianchi VIIq are spatially flat, 
Bianchi IX is positively curved and Bianchi types V and VII/^ have negative spatial 
curvature. The "open" (i.e. negatively curved) cases are of particular interest as they 
permit the focussing of anisotropic patterns into small regions of the sky. The scalar 
curvature R of the spatial sections which is given in terms of the Bianchi parameters as 
the following convenient form 

R = ^(2nin2 + 2nin3 + 2n2n3 — — nl — nl) — Qo? (11) 

For Bianchi V we have ni = n2 = na = so that R = — 6a^. In Bianchi VII/^ we have 
ni = but n2 7^ and 7^ 0; the parameter h is defined hy h = 0? /n2n-^^ it is related 
to the parameter x which defines the 'spiralness' of the temperature patterns. 

For reference, the scalar curvatures of all the models we discuss hereafter are shown in 
Table 2. 

3. Bianchi Cosmologies 

3.1. Basics 

The models we consider are based on Einstein's general theory of relativity and we use 
the field equations in the form 

Gab = Rah — -Rgah Tab ^Qab, (13) 
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with Rab being the Ricci Tensor, R the Ricci scalar, Tab the energy-momentum tensor 
and A the cosmological constant. Indices a and b run from to 3 (c.f. i, j and a, /? 
which run from 1 to 3). We use units where SirG — c — 1. In terms of a coordinate 
system x^, the metric Qab is written 

ds^ = gabdx^'dx^ = {hab - UaUb)dx''dx^ , (14) 

where is the fluid velocity; the signature of Qab is ( — h ++)• As we have already 
explained, the components of the metric Qab describing a Bianchi space are invariant 
under the isometry generated by infinitesimal translations of the Killing vector fields. 
In other words, the time-dependence of the metric is the same at all points. The Einstein 
equations relate the energy-momentum tensor Tab to the derivatives of Qab so if the metric 
is invariant under a given set of operations then so are the physical properties encoded 
in Tab. 

Before proceeding further, let us comment further on the degree of generality of the 
various Bianchi models which we touched on in the previous section. An alternative way 
to quantify this, rather than looking at the group structure, is to work out the number of 
arbitrary constants needed to specify the solutions. This seems more interesting from a 
physical point of view, as we are interested in the solutions to the field equations rather 
than the groups themselves. The number of arbitrary constants depends on the form 
of the energy-momentum tensor. In Table 3 we give the results for vacuum and perfect 
fiuid equations of state. The appearance of /i = — 1/9 as a special case in this table 
relates to the fact that two of the Einstein constraint equations become null identities for 
this particular choice of h. From the Table it emerges that the "most general" vacuum 
solutions are types VII/^, VI/^, VIII, Ylh=-i/^ and IX, all of which have four arbitrary 
parameters. The least general is the Bianchi Type I vacuum solution, which has only 
one. 

For cases describing perfect fiuids the situation is a little more subtle. One would 
expect to have four additional parameters to specify these compared to the vacuum 
solutions, but the table shows that 5 = r + 4 is not always the case. This is the case 
because the Einstein equations place additional restrictions on the form of Tab allowed 
in Types I and II. For example, if a perfect fluid is added to the vacuum Type I solution 
then the form of the metric requires all the time-space components of the Ricci tensor 
to be identically zero. This means that the energy-momentum tensor Tab must have 

To. = (15) 

for a = 1, 2 or 3. This in turn means that the matter must be comoving, i.e. its 
velocity is = 5^. Only one free parameter is therefore needed to specify the solution, 
the energy density ji. The perfect fluid case of Bianchi Vl/^^-i/g is also peculiar, in that 
it is not as general as Bianchi VI/^, VII/^, VIII or IX because the degeneracy described 
above only appears in vacuum. 
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Table 3. The Bianchi types shown in terms of the group dimension p and number of 
arbitrary constants needed to specify the model on a given constant time surface, in 
vacuum (r) and with a perfect fluid equation of state {s). 



3.2. Example: The Kasner Solution 

General solutions in closed form of the Einstein equations are only known for some 
special cases of the Bianchi types, which demonstrates the difficulty of finding 
meaningful exact solutions in situations of restricted symmetry. There is, however, 
one very well known example - the Kasner solution [74j- which is a useful illustration of 
the sort of behaviour one can obtain and which therefore provides a useful pedagogical 
route into a more general treatment of Bianchi cosmologies. 

The Kasner metric, which describes a space belonging to Bianchi Type I, has the 

form 

ds^ = dt^ - X^{t)dxl - Xldxl - Xldxl (16) 

Substituting this metric into the Einstein equations (with A = and a perfect fiuid with 
pressure p and density ji) yields 

in which = X1X2X3. Note that this emerges from the diagonal part of the Einstein 
equations so the summation convention does not apply. One also obtains 



X1X2 X2XS X^Xi 
This is easy to interpret: the spatial sections expand at a rate Xi/Xi in each direction. 
The mean rate of expansion is just 

d 1 fx, , X2 , XA 



a 3 \ Xi X2 X^ 



(19) 
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In the neighbourhood of an observer at the centre of a coordinate system Xi, fluid 
particles will move with some velocity Ui. In general, 

du^_l(du^_ dvA ,\(dih dvA _ 

dxj ~ 2 [dx, dxj^ 2 [dx, ^ dxj~ ^ ^ ' 

where uJij is the rate of rotation: in more famihar language, the vorticity vector 
= ^ijk^jk which is just the curl of Ui. The tensor Oij can be decomposed into a 
diagonal part and a trace-free part according to 

9,j = h,j9 + a,j, (21) 

where an = 0. In this description 9^ aij and uJij respectively represent the expansion, 
shear and rotation of a fluid element. In the Kasner model we have 

9 = 3{d/a) (22) 

and 

uJij = 0. (23) 

As we shall see below, more complicated Bianchi models have non-zero rotation. We 
can further write evolution equations for 

In particular we get 

+ 9a, = (25) 
which can be immediately integrated to give 

a^ = ^, (26) 

where the are constants such that Si + S2 + S3 = 0. The Kasner solution itself is for 
a vacuum p = /x = 0, which has a particularly simple behaviour described by Xi = Ait^ 
where pi +p2 +P3 = pi +p\ +p\ = 1. Notice that in general these models possess a shear 
that decreases with time. They therefore tend to behave more like an FRW model as 
time goes on. Their behaviour as t ^ is, however, quite complicated and interesting. 

3.3. Tetrad Frame 

It is convenient to follow ^Oj, introducing a tetrad basis constructed From a local 
coordinate system by 

d_ 

such that 



ea = (27) 



gab = eielgij = e\eu = diag(-l, +1, +1, +1) (28) 
meaning that the tetrad basis Ba is orthonormal. The Ricci rotation coefficients, 
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are the tetrad components of the Christoffel symbols; semicolons denote covariant 
derivatives. In general, the operators defined by equation f l27l ) do not commute: they 
generate a set of relations of the form 

[6^,65] =Yab^c- (30) 

These Ricci rotation coefficients are just 

^abc = ^{labc + 7ca6 " Ibca) - (31) 

The matter ffow is described in terms of the expansion i^ab and shear aab- 

Ua-b = ^ab + " 

C^ab = ^ab - ^"^hab = " Hhab, (32) 

where i9 = Ti^'dab) = '^aa and the magnitude of the shear is = a^^Gabl"^' We now 
take the time-like vector in our basis to be the ffuid ffow velocity so that = Sq and 
Ua = —^a- The remaining space-like vectors form an orthonormal triad, with a set 
of commutation relations like that shown in equation f l30l ) but with an explicit time 
dependence in the "structure constants" describing the spatial sections: 

[ei,e,]=7£(t)efc. (33) 

Without loss of generality we can write 

7^J = ^^Jl^^^ + ^J^^ " ^^^^J. (34) 

for some tensor n^j and some vector a^. The Jacobi identities require that Uija^ = so 
we choose = (a, 0,0) and Uij = diag(ni, 77,2, 713). The four remaining free parameters 
are used to construct the Bianchi classiffcation described brieffy above, and more in 
detail elsewhere ^ SQI [TOl [711 [72l [73] . 



4. Radiation Transport in Anisotropic Cosmologies 

Having established some general results about Bianchi models, we now turn to the 
problem of calculating the temperature and polarization patterns they produce. Our 
general method is to generalize the Liouville equation so that it can be applied to a 
complete description of photons travelling through a curved space-time. This requires 
that we set up radiation distribution functions that incorporate all the Stokes parameters 
needed to speciffc polarized radiation. We also need to include a source term that 
describes the effects of Thomson scattering by free electrons for the entire history from 
decoupling to the observed epoch. In doing this we follow closely the methods of [42j 
and also 03]. Note, however, that they use a different deffnition of the Ricci rotation 
coefficients. 
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Radiation Description 

Our expansion of the distribution functions into multipoles is based on the usual 
Stokes parameters, and on spin- weighted spherical harmonics. The transfer equation 
for polarized radiation propagating through space-time can be described by a (complex) 
photon distribution comprising components with spin-weights and 2, i.e. 

Here /, Q, U and V are the usual Stokes parameters that describe polarized radiation. 
As we shall see, however, V (which measures circular polarization) does not arise in this 
context. The degree of (linear) polarization is defined by 

P=^3^ (36) 

The polarization orientation is described by an angle v, where 

.^\.^. 

The relativistic form of the Liouville equation is 

The unpolarized part of the radiation distribution is described by spin-zero quantities, 
N (i.e. quantities invariant with respect to rotations around the ray direction k^). In 
terms of an afiine parameter A along the photon path one calculates the total change of 
N as 

dN _ ^dN dN dp"" 

and the photon path p^{\) is determined by the geodesic equation expressed in the 
tetrad notation we introduced above as 



so that 



^ = -r^„,py = (40) 

in which p^ = (e^ek^). l^ are space-time unit vectors and e is energy of photon. Where 
the 7^ is defined as 

^« = r«oo + V\k^ + T\k^ + T\k^kK (42) 

If the change in along a photon path arises from collisions only, one obtains the 
following (Boltzmann) equation 

IdN dN , ^^^p dN ^dN , y ^ .- , _. ^^^^ 

In the case of polarized radiation we need to extend the description of the radiation field 
to include both spin-0 and spin-2 components. This requires us to generalize which 
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can be decomposed into parts A^^ and A^^. In the transfer equation for A^^ these give 
rise to additional terms related to the change of angles 9 and 0, and and extra rotation 
of polarization (see the Appendix in [54j). 

^{m% + md2)N^ = 2ico89^N^ + ^{m'd + rhd)N^ (44) 
v2 dX y2 

and 

_ 2iiV2^ = -2i cot e'^N^ - 2N^m'm^e{T\^ + fc'r\). (45) 
d\ d\ 

(see [43] for details). The first two extra terms cancel out in the Liouville equation so 
we obtain the following simplified form 

+ 5\iN^{V^^^e'^^k^ + r^fc^fc^e^^^). (46) 

We used the antisymmetry of a and c in F^^ and the relation, \{m^m^ —fh}m^) = —^e^^^k^ 
to obtain the imaginary terms. The operator Va on the left hand side preserves spin 
weight {A = 0, 2); the corresponding angular operator is i9 = ^ / V2{m'^B + m^S). 

4-2. Scattering 

Photon scattering can be described by the addition of a source term to the right hand 
side of the Liouville equation so that it can be written 

VaN^ = t{-N^ + Ja), (47) 

in which Ja takes the form 

r - dVt' 

Ja = J [pABie, 0, 9', m^^io'. 4>') + pAB{e, 0, e\ <p')n''{9', <j>')]— (48) 

in terms of the scattering matrices pab and Pab- 

Poo = ^ + lk"'k"'' + Ik'k'' 
Z o Z 

Po2 = Po2 = -^k m'" 

Pm = \ + \k'''k^^' - \k'k' 
3 / 

P22 = -m'^m'^\ (49) 

The emission term J a contains only harmonics up to / = 2, since all other terms vanish 
in virtue of the orthogonality relations for spherical harmonics 

Jo = ReiV° + IknmN^ + ^k^'ReN^ - ^k^'ReNl (50) 
J,= - im^'^ReiVO + ^mHeNl (51) 
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The radiation modes with / < 2 are damped as weU as re-radiated by Thomson 
scattering, while higher-order modes / > 2 are only damped. 

4 . 3. Transfer Equation 

We now expand the distribution function in our Boltzmann equation in terms of 
multipole components: 

A^o = + N^k' + N^jk'^ + . . . , 

= Nf.m'^ + . . . (52) 
where A;* expresses the three-dimensional ray direction, i.e. 

k'^ = (cos 6, sin 6 cos 0, sin 6 sin 0) 

j^ik _ j^ij^k _ l^ifc 
3 

1 ,dk^ i dk\ 
m — ( \ • I 

^^de sin^ d(t)^ 

m'^ = rn'm!" (53) 

in which rrfki — 0, m}fhi — k'^ki = 1 and the bar indicates complex conjugate. The 
number of indices of the k and m polynomials characterizes the multipole order of the 
corresponding contributions to anisotropy and polarization. This leads to the following 
equations for the moments of these quantities 

47r 3 

47r 15^ 3 ^ 

47r 10^ ^ 15 

47r 30^ V / 

where the integration is taken over the unit sphere. 

We are now in a position to write down equations for the evolution of the 
components of the distribution function, as follows: 



I- 



I 
I 



^*A^S - l^kd.Nf, + k^^Nl = -^(^ReiVS + ^ReiVg + zlmiVj). (55) 

We have introduced the differential operators dt — d/dt and — ed/de. The latter is 
used in the definition of C,\ 

a(lniVO) _^^ a(lniV°) 
9 In £ de 
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The coefficients arising in the transfer equation can be obtained in the foUowing form 

"^kl k k k 

Ci = - - (r + r ^^6ii + 3r qqSu) 

"fc 

Aj = 3^ QiSij - -T Q-Sij - -T QjSii 
- 3^ ii^ij - y ij - y ji 

^ kl l 

Fij = j^iSij - -{SkiSijT + r j^-Sji + r j^jSu) 

5 2 

- -(r ^^(^j/ + r jj^Sii) + -r rnm^ik^jl 

(57) 

Because fc^-^ and m^^ are both symmetric and traceless, it follows that N^j and N^j are 
too. The coefficients listed above must therefore be made symmetric and traceless on 
the index pairs ij and kl. 

Equations f l55l ) describe the behavior of the lowest-order angular modes of a general 
radiation field in a general space-time, with two main assumptions. First, in the 
orthonormal frame, we assume that the "spatial" derivatives ^^A^ vanish for all quantities 
of interest because of homogeneity. Second, we assume that the radiation field is 
described by a Planck distribution at all times. These assumptions reduce the set 
of equations needed to a linear system of coupled ordinary differential equations with 
"time" -dependent coefficients. The patterns that are produced therefore depend both 
on the background cosmology and the initial conditions. 

Thomson scattering does not affect the component ReN^ but it does the term 
ImA^Q that describes circular polarization. The mode describing linear polarization, Nfj^ 
is coupled to higher-order modes of the radiation field but not directly through Nq. 
Any non-zero term (7f ^ would produce a dipole variation of the radiation distribution 
represented by but no dipole can be produced this way in the particular case of 
Bianchi I (nor indeed for the FRW case). The presence of a dipole component is 
inevitable in other Bianchi models since, even with A^^^^ = initially, A^^^^ becomes 
different from zero if F^q^ 7^ or if N^iT^i- + N^j^T^n differs from zero. In a similar 
manner a quadrupole component can always be generated from the isotropic mode A'q , 
if Eik is different from zero which means that the ffuid ffow possesses some kind of 
shearing motion. 
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It is clear from the system of equations (l55l) that a gravitational field alone is not 
able to generate polarization. Initially unpolarized radiation collisionlessly propagating 
in an arbitrary gravitational field must remain unpolarized. However, a quadrupole 
mode of unpolarized radiation generates a linear polarization component at / = 2, if 
Thomson scattering is present. This is the standard mechanism by which polarization 
is generated from radiation anisotropics in the early Universe. Since a quadrupole mode 
of unpolarized radiation is generated from the isotropic component if Eik 0, the 
cosmological gravitational field could therefore be indirectly responsible for generating 
polarization if it first generates a quadrupole anisotropy. Equations f l55l ) also show that 
an initial monopole Nq produces a non-zero quadrupole A^^^- via the shear Eij = —cjij 
and A^^^^ is subsequently coupled with Nf- by C^K The effect of shear on the radiation 
is to generate a quadrupole anisotropy by redshifting the it anisotropically; in some 
models, dipole and higher order multipoles would also arise. 

To summarize, then. In order to get interesting higher-order patterns in the 
radiation background we must either have non-zero shear if there is no initial quadrupole 
or have non-zero initial quadrupole and monopole terms if there is no shear. 



4.4. Geodesic Equations 

We now follow the follow the convention in ref. [42j to construct the equations describing 
geodesies in the models we consider. From each components of eq. ( l40l) time variation 
of e and direction vector have the form such as 

^ = ^ = -eW (58) 

c^A "c^A^ ^'d\ ~ ^^^^ 
Using the relations 

dp' = k'de + edk' = k'ds + £{a'de + sin Ob'dc^) (60) 

and orthogonality of , and fc^, we obtain the following equations for the time 
variation for 9 and 0: 

d9 1 df . 

-7T = -^i-Tv = -^^a (61) 

aA e dX 

# ^ l^df ^ _^ J^y .g2) 

dX £ sin 9 dX sin 9 

since the 7^ are also represented by Ricci Coefficients F^^. Here, we also neglect the 
effect shear which would make (F^^- = 0). Finally we obtain variation terms: 

1 d9 

--— = sin ^ [a + cos sin 0(77.3 — 77,2)] (63) 
£ dX 

--^ = cos ^[ni - 7^3 + {ns - 7^2) cos^ 0] (64) 
£ dX 

The change of polarization angle ij) is expressed by 

= im'fh^Tlk^ - cot 9ba^ = -nt - nlk'^ + cos 9-^. (65) 
£dX 6 ^ " £dX ^ ^ 
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Although the Tf- contain the Bianchi vector a^, only the components appear in the 
change of since aik^k^e^^^ = by symmetry of the / and m indices. Consiquently this 
term does not give different results for the models such as type I and V or VII/^ and 
VIIq. Applying this to each Bianchi types we obtain, for Types I and V, 

For Types VII/^ and VIIq we get 

1^ = 2n(l - sin^ 6) + - cos 9^ = -ncos^ 9 (67) 
s dX £ dX 

and for Type IX we have 

Sdx = ^''^ 

Note that we consider the FRW limits i.e. Bianchi vectors as = n. The term ^ 
is important since it gives rise to mixing terms between the E and B modes in the 
Boltzmann equation. 



5. Temperature and Polarization Patterns 

In this section we will present representative examples of the temperature and 
polarization patterns produced in the models we have discussed, computed by 
numerically integrating the system of equations derived in Section 4. The patterns 
produced depend on the parameters chosen for the model in question and also, as we 
have explained in the previous section, on the choice of initial data. In the following, 
primarily pedagogical, discussion we do not attempt to normalize the models to fit 
current cosmological observations but restrict ourselves to phenomenological aspects of 
the patterns produced. The overall level of temperature anisotropy depends on the 
choice of parameters in the Bianchi models that express the extent of its departure from 
the FRW form. Since we do not tune this to observations, the amplitude is arbitrary, as 
in the cosmic epoch attributed to each of the evolutionary stages. Moreover, the overall 
degree of polarization depends strongly on the ionization history through the optical 
depth T which appears in Eq. f l47l) . We shall not attempt to model this in detail in this 
paper either. What is important, however, is that the geometrical relationship between 
the temperature and polarization patterns does not depend on these factors; it is fixed 
by the geometric structure of the model, not on its normalization. 

Of course we compute only the coherent part of the radiation field that arises from 
the geometry of the model. Any realistic cosmological model (i.e. one that produces 
galaxies and large-scale structure) must have density inhomogeneities too. Assuming 
these are of stochastic origin they would add incoherent perturbations on top of the 
coherent ones produced by the background model. 

For illustrative purposes we have chosen cases where the initial conditions produce 
a pure quadrupole anisotropy of tesseral form, i.e. and / = 2 spherical harmonic mode 
with m = 1 as shown in Figure 1. We gratefully acknowledge the use of the Healpix 
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Figure 2. Polarization maps, i.e. polarization amplitude (left), Stokes parameter Q 
(middle) and U (right) of Bianchi type I. 

software [75] in creating this and all the other all-sky maps shown hereafter. 

Other choices are, of course, possible. A quadrupole with m = would produce 
a zonal pattern, and one with m = 2 would be a sectoral mode [38j; these choices 
are discussed at some length in ref. [57]. One could also generate more complicated 
patterns by having an observer who is not at rest in the frame we are using, which would 
introduce an additional dipole anisotropy. We will not discuss this possibility further in 
this paper. 

The simplest case is obviously that of Bianchi I, but this is nevertheless of some 
interest because a Universe of this type could in principle account for the presence of 
a low quadrupole [76[ |77j . In this example the temperature pattern does not evolve at 
all with time, so one can simply treat the initial quadrupole as a free parameter. The 
polarization patterns arising in this model, which do not evolve with time either, are 
shown in Figure 2. 

We next turn our attention to Bianchi types V, Vllh and VIIq. These models have 
a single preferred axis of symmetry. The alignment of the shear eigenvectors relative to 
this preferred axis determines not only the dynamical evolution of the model through 
the field equations, but also the temperature and polarization pattern which, as we 
shall see, gets imprinted into the cosmic background radiation. Figure 3 shows (from 
bottom to top) the time evolution of the temperature pattern in these models. Note 
that, in Type V (left), the initial quadrupole retains its shape but gets focussed into a 
patch of decreasing size as time goes on. This is due to the effect of negative spatial 
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Figure 3. The time evolution of the temperature maps for Bianchi types V, Vllh and 
VIIq. Time increases from bottom to top. 

curvature. In Bianchi Type VIIq (right) the effect of rotation and shear is to twist the 
initial quadupole into a spiral shape that winds up increasingly as the system evolves. 
In the middle case, Bianchi Type VII/^ we have a combination of the two cases either 
side: there is both a focussing and a twist. This case produces the most complicated 
temperature pattern. 

In the following three figures we examine the polarization pattern produced in the 
models shown in Figure 2. First, in Figure 3, we have Bianchi V. These results show 
that while the polarization pattern alters with time in this case, its general orientation 
on the sky does not (as is the case with the temperature pattern). The implications of 
this for the production of cosmological B-mode polarization was discussed by ref. [55j. 

Figure 5 shows an example of Bianchi Type VII/^. Note the prodigious twisting of 
the polarized component of the radiation field, as well as a concentration of the degree 
of polarization defined by Equation f l36l) . This case generically produces a large amount 
of mixing between Q and U during its time evolution. 

Figure 6 shows, not unexpectedly, that Bianchi Type VIIq produces a similar 
interweaving of the Q and U configurations, but without the focussing effect in the 
total polarized fraction. 

Finally, for completeness we show results for Bianchi Type IX in Figure 7. 
This provides an interesting example of polarization behaviour because, in terms of 
polarization degree, the pattern does change at all with time but the Stokes maps Q 
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Figure 4. The time evolution of polarization maps for Bianchi V: degree of linear 
polarization (left); Stokes parameter Q (middle); and U(riglit). Time increases from 
bottom to top. 




Figure 5. Same as Figure HI but for Bianchi type Vllh- 
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Figure 6. Same as Figure HI but for Bianchi type VIIq. 



and IJ do evolve. Bianchi Type IX models have positively curved spatial sections, and 
are equivalent in some sense to FRW models with the addition of circularly-polarized 
gravitational waves. These cause a rotation of the polarization angle but do not change 
the overall magnitude. The positive curvature does not allow for any focussing effects 
so the temperature pattern does not change with time either. 



6. Discussion and Conclusions 



The goal of this work was to compute the temperature and polarization patterns 
produced in anisotropic relativistic cosmologies described by various Bianchi types. 
We focussed on those types that contain the standard homogeneous and isotropic 
FRW background as a limiting case. We constructed an appropriate description of the 
radiation field in terms of spin-0 and spin-2 components representing the unpolarized 
and polarized parts, respectively. We integrated the equations for the radiation field 
numerically, and presented illustrative examples here. 

The basic point that emerges from this study is related to the physical origin 
of CMB polarization: radiation affected by Thomson scattering from an electron 
in a radiation field possessing a quadrupole anisotropy will inevitably be partially 
linearly polarized. In the context of standard cosmological models, the environment 
of different electrons varies owing to the presence of density inhomogeneities and a 
background of gravitational waves. These sources of variation are stochastic so the 
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Figure 7. Same as Figure IH but for Bianchi type IX. 

variations in the polarized component of the radiation field, though correlated with 
the temperature variations, are essentially incoherent. In Bianchi cosmologies, however, 
global homogeneity requires that each electron sees the same quadrupole. The polarized 
part of the radiation field is therefore coherent, and is in a fixed relationship to the 
temperature variation (once the model is specified). 

We have exploited the formalism used to generate the examples shown here to 
produce an extensive test-bed of CMB temperature and polarization maps derived from 
Bianchi universe models of the types discussed in this paper. In future work we will 
analyse these maps using a variety of statistical measures of anisotropy. It is important 
to understand how standard statistical techniques, designed to be applied to stationary 
stochastic fluctuations, perform when applied to patterns which are neither stochastic 
nor stationary. Our Bianchi models, for example, yield patterns which have highly 
correlated spherical harmonic components in contrast to the case of a homogeneous 
Gaussian random field in which the harmonic modes are independent. It is important 
to understand how CMB fluctuations arising from Bianchi-type universes impact on 
orthodox analysis procedures and whether they produce characteristic signatures when 
analysed in this way, particularly in the presence of additional stochastic fluctuations. 
Higher-order statistics will be necessary to provide a full characterization of the coherent 
anisotropic fluctuations we have studied here. Moreover, present and future CMB 
experiments are largely aimed at improving the precision of polarization measurements 
in order to find evidence of a stochastic background of primordial gravitational waves. 



Temperature and Polarization Patterns in Anisotropic Cosmologies 



22 



Techniques advocated to measure polarization pattern over the celestial sphere are 
therefore also generally tuned to detect signals of an incoherent nature. The presence 
of coherent signals in the CMB could be an indication of physics beyond the standard 
cosmological framework. This will hopefully lead us to better ways of searching for 
departures from the concordance cosmology using the next generation of datasets. 
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